Lecture 15: Learning probabilistic models 

Roger Grosse and Nitish Srivastava 


1 Overview 

In the first half of the course, we introduced backpropagation, a technique we used to 
train neural nets to minimize a variety of cost functions. One of the cost functions we 
discussed was cross-entropy, which encourages the network to learn to predict a probability 
distribution over the targets. This was our first glimpse into probabilistic modeling. But 
probabilistic modeling is so important that we’re going to spend almost the whole second 
half of the course on it. This lecture introduces some of the key principles. 

Actually, there’s aren’t any major new ideas in this lecture. You’ve already seen maxi¬ 
mum likelihood estimation in the context of neural probabilistic language models (Coursera 
Lecture D, in-class Lecture 7, and Assignment 1). Coursera Lecture J introduced the full 
Bayesian approach and the maximum a-posteriori (MAP) approximation. All we’re doing 
here is stating these principles in slightly more general terms, and working through lots of 
examples in order to gain a better intuition. Once you’ve gotten more practice with these 
techniques, it’s a good idea to go back and revisit those lectures. 

This lecture and the next one aren’t about neural nets. Instead, they’ll introduce the 
principles of probabilistic modeling in as simple a setting as possible. Then, starting next 
week, we’re going to apply these principles in the context of neural nets, and this will result 
in some very powerful models. 

1.1 Learning goals 

• Know some terminology for probabilistic models: likelihood, prior distribution, poste¬ 
rior distribution, posterior predictive distribution, i.i.d. assumption, sufficient statis¬ 
tics, conjugate prior 

• Be able to learn the parameters of a probabilistic model using maximum likelihood, 
the full Bayesian method, and the maximum a-posteriori approximation. 

• Understand how these methods are related to each other. Understand why they tend 
to agree in the large data regime, but can often make very different predictions in 
the small data regime. 


1 



2 Maximum likelihood 


The first method we’ll cover for fitting probabilistic models is maximum likelihood. In 
addition to being a useful method in its own right, it will also be a stepping stone towards 
Bayesian modeling. Actually, you’ve already done maximum likelihood learning in the 
context of the language model from Assignment 1. All we’re doing now is presenting the 
more general framework. 

Let’s begin with a simple example: we have flipped a particular coin 100 times, and it 
landed heads Njj = 55 times and tails Nt = 45 times. We want to know the probability 
that it will come up heads if we flip it again. We formulate the probabilistic model: 

The behavior of the coin is summarized with a parameter 9, the probability 
that a flip lands heads (H). The flips T> = - • ■,x^ 100 )) are independent 

Bernoulli random variables with parameter 6. 

(In general, we will use T> as a shorthand for all the observed data.) We say that the indi¬ 
vidual flips are independent and identically distributed (i.i.d.); they are independent 
because one outcome does not influence any of the other outcomes, and they are identically 
distributed because they all follow the same distribution (i.e. a Bernoulli distribution with 
parameter 9). 

We now define the likelihood function L{6 ), which is the probability of the observed 
data, as a function of 9. In the coin example, the likelihood is the probability of the 
particular sequence of H’s and T’s being generated: 

L{9) = p(V) = 9 Nh ( 1 - 9) Nt . 

Note that L is a function of the model parameters (in this case, 9), not the observed data. 

This likelihood function will generally take on extremely small values; for instance, 
L(0.5) = 0.5 100 ~ 7.9 x 10 -31 . Therefore, in practice we almost always work with the 

log-likelihood function, 

£(d) = log L{9) = N h log 9 + Nt log(l — 9). 

For our coin example, £(0.5) = logO.5 100 = 100 log 0.5 = —69.31. This is a much easier 
value to work with. 

In general, we would expect good choices of 9 to assign high likelihood to the observed 
data. This suggests the maximum likelihood criterion: choose the parameter 9 which 
maximizes i{9). If we’re lucky, we can do this analytically by computing the derivative and 
setting it to zero. (More precisely, we find critical points by setting the derivative to zero. 
We check which of the critical points, or boundary points, has the largest value.) Let’s try 
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this for the coin example: 


^ ^ (N h log 0 + N t log(l - 6)) 

Nn Nt 

= ~e~ ~ 1 - e 


Setting this to zero, we find the maximum likelihood estimate 


#ML 


N h 

N h + N t ’ 


i.e. the maximum likelihood estimate is simply the fraction of flips that came up heads. 
(We put a hat over the parameter to emphasize that it’s an estimate.) Hopefully this seems 
like a sensible guess for 6. Now let’s look at some more examples. 


Example 1. Suppose we are interested in modeling the distribution of temperatures 
in Toronto in March. Here are the high temperatures, in Celsius, from the first week 
of March 2014: 


-2.5 -9.9 -12.1 -8.9 -6.0 -4.8 2.4 


Call these observations x^\ ... ,x^ N \ where N = 7. In order to formulate a prob¬ 
abilistic model, we first choose a parametric form for the distribution over temper¬ 
atures. Often we choose a Gaussian distribution, not because we believe it’s an 
especially good model, but because it makes the computations easy. So let’s assume 
the temperatures are drawn from a Gaussian distribution with unknown mean and 
known standard deviation a = 5. Our likelihood function is given by: 



i=l 


Since // can take any possible real value, the maximum must occur at a critical point, 
so let’s look for critical points. Setting the derivative to 0, 

i N i 
>rr 2 


o = ' V T (l (i) _ „)2 

dn 2a 2 du v 
1=1 
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Therefore, J2iLi x ^ — p = 0, and solving for p, we get p = The 

maximum likelihood estimate of the mean of a normal distribution is simply the 
mean of the observed values, or the empirical mean. Plugging in our temperature 
data, we get /1 ml = —5.97. 

Example 2. In the last example, we pulled the standard deviation a = 5 out of a 
hat. Really we’d like to learn it from data as well. Let’s add it as a parameter to the 
model. The likelihood function is the same as before, except now it’s a function of 
both p and a, rather than just p. To maximize a function of two variables, we find 
critical points by setting the partial derivatives to 0. In this case, 


dp a 2 


i= 1 

n - — - — ' N 
<j “ da - 




N 


-\ log 2 tt - log a - ^ (x W - pf 

(x® - pf 


U=l 


da 
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da 2 a 


a a- J 
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2=1 


From the first equality, we find that pml = j y HZiLi x ^ is the empirical mean, just 

as before. From the second inequality, we hnd cjml = \J jj — A 4 ) 2 - other 

words, a ml is simply the empirical standard deviation. In the case of the Toronto 
temperatures, we get pml = —5.97 (as before) and <tml = 4.55. 

Example 3. We’ve just seen two examples where we could obtain the exact max¬ 
imum likelihood solution analytically. Unfortunately, this situation is the exception 
rather than the rule. Let’s consider how to compute the maximum likelihood estimate 
of the parameters of the gamma distribution, whose PDF is defined as 

p(x) = 

F(o) 

where T(a) is the gamma function, which is a generalization of the factorial function 
to continuous values. 1 The model parameters are a and b, both of which must take 

lr The definition is T(t) = J 0 °° x t ~ 1 e~ x da;, but we’re never going to use the definition in this class. 
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positive values. The log-likelihood, therefore, is 

N 

£(a , b) = ^ alogb — logT(a) + (a — 1) logx*-^ — bx ^ 

1=1 

N N 

= Na log b — N logT(a) + (a — 1) ^ logx^ — 6^ x^. 

i=1 i =1 

Most scientific computing environments provide a function which computes log T(a). 
In SciPy, for instance, it is scipy. special. gammaln. 

To maximize the log-likelihood, we’re going to use gradient ascent, which is just like 
gradient descent, except we move uphill rather than downhill. To derive the update 
rules, we need the partial derivatives: 

d N 

N log b — N — log T (a) + ^ log x ^ 

i= 1 
N 

iv“ wo. 

b i t 

Our implementation of gradient ascent, therefore, consists of computing these deriva¬ 
tives, and then updating a <— a + a 9^ and b <— b + where a is the learning rate. 
Most scientific computing environments provide a function to compute ^ log T(a); 
for instance, it is scipy. special. digamma in SciPy. 

Here are some observations about these examples: 

• In each of these examples, the log-likelihood function l decomposed as a sum of terms, 
one for each training example. This results from our independence assumption. Be¬ 
cause different observations are independent, the likelihood decomposes as a product 
over training examples, so the log-likelihood decomposes as a sum. 

• The derivatives worked out nicely because we were dealing with log-likelihoods. Try 
taking derivatives of the likelihood function L(6), and you’ll see that they’re much 
messier. 

• All of the log-likelihood functions we looked at wound up being expressible in terms 

of certain sufficient statistics of the dataset, such as X^i[ x ^] 2 , or 

log x'W. When we’re fitting the maximum likelihood solution, we can forget the 
data itself and just remember the sufficient statistics. This doesn’t happen for all 
of our models; for instance, it didn’t happen when we fit the neural language model 


cM 
da 

d£ 

db 
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in Assignment 1. However, it does happen for many of the distributions commonly 
used in practice. 2 

• We made a lot of questionable assumptions in formulating these models. E.g., we 
assumed that temperatures on different days were independent; in practice, the tem¬ 
perature tomorrow is likely to be similar to the temperature today. This is also true 
of models we fit previously; e.g., the Markov assumption we used to justify our neural 
language model is clearly bogus. If this were a statistics class, we’d talk about ways 
to test your modeling assumptions. But because this is a machine learning class, 
we’ll throw caution to the wind and fit models that we know are wrong. Hopefully 
they’ll still be good enough that they can make sensible predictions (in the supervised 
setting) or reveal interesting structure (in the unsupervised setting). 

2.1 Beware of data sparsity 

Maximum likelihood is a very powerful technique, but it has a pitfall: if you have too 
little training data, it can seriously overfit. The most severe pathology is when it assigns 
probability 0 to things that were never seen in the training set, but which still might 
actually happen. For instance, suppose we flip a coin twice, and it lands H both times. 
The maximum likelihood estimate of 9, the probability of H, would be 1. But this is pretty 
extreme — effectively we’re considering it impossible that the coin will ever come up T! 
This problem is known as data sparsity. 

This problem isn’t so different in principle from examples of overfitting which we dis¬ 
cussed for other loss functions. We would like our model to generalize well to data it 
hasn’t seen before, i.e. assign the new data high likelihood. We can measure the general¬ 
ization performance by holding out a separate test set, and measuring the log-likelihood 
on this test set at the very end. (As before, if we want to choose model hyperparamters, 
we’d hold out a separate validation set.) In our coin example, if we choose 9 = 1 and the 
coin comes up T even a single time in the test set, this will give us a test log-likelihood of 
—oo. Clearly it’s a bad idea to assign any outcome probability 0 if it might ever occur. 

Last week, we talked about regularization as a way to attenuate overfitting. Would that 
work here? One naive approach would be to add an L -2 penalty, —\9 2 , to the objective 
function. (We subtract the penalty since we’re maximizing.) But this isn’t quite what 
we want: it would allow (in fact, encourage) the degenerate solution 9 = 0. Instead, let’s 
look at Bayesian techniques for parameter estimation. These techniques will turn out to 
be closely related to regularization. 

2 If you’re interested in learning more, the families of distributions whose likelihoods can be written in 
terms of sufficient statistics are known as exponential families. 
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3 Bayesian parameter estimation 


In the maximum likelihood approach, the observations (i.e. the xf s) were treated as random 
variables, but the model parameters were not. In the Bayesian approach, we treat the 
parameters as random variables as well. We define a model for the joint distribution 
p(9,V) over parameters 9 and data V. (In our coin example, 9 would be the probability 
of H, and V would be the sequence of 100 flips that we observed.) Then we can perform 
the usual operations on this joint distribution, such as marginalization and conditioning. 

In order to define this joint distribution, we need two things: 

• A distribution p(9), known as the prior distribution. It’s called the prior because 
it’s supposed to encode your “prior beliefs,” i.e. everything you believed about the 
parameters before looking at the data. In practice, we normally choose priors to be 
computationally convenient, rather than based on any sort of statistical principle. 
More on this later. 


• The likelihood p(V \ 9), the probability of the observations given the parameters, 
just like in maximum likelihood. 


Bayesians are primarily interested in computing two things: 


• The posterior distribution p(9\V). This corresponds to our beliefs about the 
parameters after observing the data. In general, the posterior distribution can be 
computed using Bayes’ Rule: 


p(9\V) 


p(9)p(V | 0) 
f p(9')p(V | O') cW ' 


However, we don’t normally compute the denominator directly. Instead we work with 
unnormalized distributions as long as possible, and normalize only when we need to. 
Bayes’ Rule can therefore be written in a more succinct form, using the symbol oc to 
denote “proportional to”: 

p(9 | V ) oc p(9)p(V | 9). 


• The posterior predictive distribution p(T>' \T>), which is the distribution over 
future observables given past observations. For instance, given that we’ve observed 
55 H’s and 45 T’s, what’s the probability that the next flip will land H? We can 
compute the posterior predictive distribution by computing the posterior over 9 and 
then marginalizing out 9: 


Pip | V) 


J p(9 | V)p{V' | 9) d 9. 
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Figure 1: The PDF of the beta distribution for various values of the parameters a and b. 
Observe that the distribution becomes more peaked as a and b become large, and the peak 
is near a/(a + b). 


3.1 The full Bayesian approach 


Let’s figure out the posterior distribution and posterior predictive distribution for our 
coin example. We’ve already specified the likelihood, so it remains to specify the prior. 
One option is to use an uninformative prior, which assumes as little as possible about 
the problem. In the case of the coin, this might correspond to the uniform distribution 
p(9) = 1. (There is no single recipe for choosing an uninformative prior; statisticians have 
a few different recipes which often, but not always, agree with each other.) 

Alternatively, we can draw upon our lifetime of experience flipping coins. Most coins 
tend to be fair, i.e. the come up heads around 50% of the time. So perhaps our prior 
should make 9 = 0.5 more likely. There are a lot of distributions which can do this, but a 
particularly useful one is the beta distribution, parameterized by a,b > 0, and defined 
as: 


p{6-,a,b) 


r(a + b) 
r(a)r(6) 


9 a ~ x 


(1 -Of- 1 . 


This distribution is visualized in Figure 1. Why did we choose the beta distribution, of all 
things? Once we work through the example, we’ll see that it’s actually pretty convenient. 
Observe that the first term (with all the T’s) is just a normalizing constant, so it doesn’t 
depend on 6. In most of our computations, we’ll only need to work with unnormalized 
distributions (i.e. ones which don’t necessarily integrate to 1), so we can drop the ugly 
normalizing constant and write 


p{9-,a,b) ^9 a - 1 {l^9) b - 1 . 













A few values are plotted in Figure 1. From these plots, we observe a few things: 

• The distribution appears to be centered around a/(a + b). In fact, it’s possible to 
show that if 6 ~ Beta(a, b), then E[0] = a/(a + b). 

• The distribution becomes more peaked for larger values of a and b. 

• The values a = b = 1 correspond to the uniform distribution. Therefore, we can 
simply treat the uniform prior as a special case of the beta prior. 


Now let’s compute the posterior and posterior predictive distributions. When we plug 
in our prior and likelihood terms for the coin example, we get: 


P {e | v) oc p(o)p(v | e) 

\e a - 1 (i-d ) b ~ 1 


oc 


0 Nh (1 - 0 ) Nt 


_ Qa-1+N H ^ _ Q^b—1+N T 


(1) 

( 2 ) 

(3) 


But this is just a beta distribution with parameters Nh + a and Nt + b. Let’s stop and 
check if our answer makes sense. As we observe more flips, Nh and Nt both get larger, 
and the distribution becomes more peaked around a particular value. Furthermore, the 
peak of the distribution will be near Nh/(Nh + Nt), our maximum likelihood solution. 
This reflects the fact that the more data we observe, the less uncertainty there is about 
the parameter, and the more the likelihood comes to dominate. We say that the data 
overwhelm the prior. 

We now compute the posterior predictive distribution over the next flip x 

0pred = Pl V = H \V) 

= J p{6 | V)Pr(x' = H \ 0)6.6 

= J Beta(0; N H + a, N T + b) ■ 6 d 9 

= ®Beta(6»;Af /i -+a,Af T +6)[^] 

_ N h + a 
Nh + Nt + a + b ’ 

where the last line follows from the formula for the expectation of a beta random variable. 
Again, let’s check if this is reasonable. If Nh and Nt are large, this is approximately 
Nh/(Nh+Nt), so our predictions are similar to the ones we get using maximum likelihood. 
However, if Nh and Nt are relatively small, then the predictive distribution is smoothed, 
i.e. less extreme than the maximum likelihood one. The value # pre d is somewhere in between 
the prior and the maximum likelihood estimate. 
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Figure 2: Plots of the prior, likelihood, and posterior for the coin flip example, with the 
prior Beta(2,2). (Left) Small data setting, Nh = 2, Nt = 0. (Right) Large data 
setting, Nfj = 55, Nt = 45. In this case, the data overwhelm the prior, so the posterior is 
determined by the likelihood. Note: for visualization purposes, the likelihood function is 
normalized to integrate to 1, since otherwise it would be too small to see. 

OK, back to an earlier question. Where did our choice of prior come from? The key 
thing to notice is Eqn 3, where the posterior wound up belonging to the same family of 
distributions as the prior. Why did this happen? Let’s compare the formulas for the beta 
distribution and the likelihood: 

p(9) = Beta(0; a, b ) oc Q a ~ l ( 1 — 0) 6_1 
p(V\6) oc 6 Nh ( 1 - 9) Nt 

In other words, the prior was chosen to have the same functional form as the likelihood. 3 
Since we multiply these expressions together to get the (unnormalized) posterior, the pos¬ 
terior will also have this functional form. A prior chosen in this way is called a conjugate 
prior. In this case, the parameters of the prior distribution simply got added to the 
observed counts, so they are sometimes referred to as pseudo-counts. 

Let’s look at some more examples. 

Example 4. Let’s return to our problem of estimating the mean temperature in 
Toronto, where our model assumes a Gaussian with unknown mean p and known 
standard deviation a = 5. The first task is to choose a conjugate prior. In order to 
do this, let’s look at the PMF of a single data point: 

( | \ 1 ( {x-p) 2 \ 

3 The oc notation obscures the fact that the normalizing constants in these two expressions may be 
completely different, since p(6) is a distribution over parameters, while p(V \ 6) is a distribution over observed 
data. In this example, the latter normalizing constant happens to be 1, but that won’t always be the case. 
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If we look at this as a function of ji (rather than x), we see that it’s still a Gaus¬ 
sian! This should lead us to conjecture that the conjugate prior for a Gaussian is a 
Gaussian. Let’s try it and see if it works. 

Our prior distribution will be a Gaussian distribution with mean p, pr i and standard 
deviation cr pr i. The posterior is then given by: 
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The last step uses a technique called completing the square. You’ve probably done 
this before in a probability theory class. So the posterior distribution is a Gaussian 
with mean /i pos t and standard deviation ( 7 post . 

The formulas are rather complicated, so let’s break them apart. First look how (J post 
changes if we vary the prior or the data. 


— As we increase the number of observations N, the denominator gets larger, so 
(jpost gets smaller. This should be intuitive: as we observe more data, the pos¬ 
terior gets more peaked, which corresponds to the posterior standard deviation 
decreasing. 
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— What if we increase the prior standard deviation cr pr i or the observation standard 
deviation cr? Then the denominator gets smaller, which means cr post gets larger. 
This should be intuitive, because increasing the uncertainty in either the prior 
or the likelihood should increase the uncertainty in the posterior. 

Now let’s look at the formula for p pos t. It takes the form of a weighted average 
of the prior mean p pr ; and the maximum likelihood mean / T By weighted 

average, I mean something of the form 

ar + b-s 
a + b 

where the weights a and b are both positive. This is a weighted average of r and s; if 
a is larger, it is closer to while if b is larger, it is closer to s. For p pos t, the weights 
for the prior mean and maximum likelihood mean are l/<7p ri and N/a 2 , respectively. 
Let’s see what happens if we change the problem. 

— As the number of observations N gets larger, the weighted average becomes 
closer and closer to the maximum likelihood estimate. This should make sense: 
as we get more information, our prior beliefs become less relevant. 

— If we increase the prior standard deviation cr pri , then 1/crjG gets smaller, so 
there is less weight on the prior. On the other hand, if we increase the standard 
deviation cr of the observations, then N/a 2 gets smaller, and there is less weight 
on the maximum likelihood solution. In other words, whichever of the two terms 
we are more confident about should count for more. 

Observe that YliLi x ^ is a sufficient statistic, since it is the only thing we need to 
remember about the data in order to compute the posterior. 

Finally, let’s take a look at the posterior predictive distribution. We compute this as 
p(x' | V) = J p(p | T>)p(x' | p) dp 

= J Gaussian(/Lt; p pos t, <J P ost)Gaussian(a; / ; p, a) dp 
= Gaussian(x / ; ^ post , y^post + cr 2 ) 

The last step uses the formula for the convolution of two Gaussian distributions. Now 
let’s see how it behaves. 

— When there are no observations (i.e. N = 0), then p pos t and <7 post are the prior- 
mean and standard deviation. The predictive distribution is centered at p pos t, 
but more spread out than the prior. 
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Figure 3: The prior, posterior, and posterior predictive distributions for the Toronto tem¬ 
peratures example. 


— When N is very large, the mean of the predictive distribution is close to the 
maximum likelihood mean, and the standard deviation is very close to u. In 
other words, it makes almost the same predictions as the maximum likelihood 
estimate. 

The prior, posterior, and posterior predictive distributions are all shown in Figure 3. 

For both the coin and Gaussian examples, the posterior predictive distribution had 
the same parametric form as the model. (I.e., it was a Bernoulli distribution for the coin 
model, and a Gaussian distribution for the Gaussian model.) This does not happen in 
general; often the posterior predictive distribution doesn’t have a convenient form, which 
is part of what makes the full Bayesian approach difficult to apply. 

3.2 The difficulty of the full Bayesian approach 

We’ve seen two different ways to learn the parameters of a probabilistic model. Maximum 
likelihood is based on optimization, while the full Bayesian approach is based on computing 
integrals. In either case, for some of the commonly used distributions, we can derive a 
closed-form solution. However, for many important models (such as multilayer neural 
nets), there’s no closed-form solution. As we saw in Example 3, if we can’t find a closed 
form, we can still maximize the log-likelihood using gradient ascent. 

But for the Bayesian approach, we need to compute an integral in order to marginalize 
out the model parameters. If we only have a few parameters, we can do this using nu¬ 
merical quadrature methods. Unfortunately, these methods are exponential in the number 
of variables being integrated out. If we’re trying to fit a neural net with thousands (or 
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even millions) of parameters, this is completely impractical. There are other methods for 
integration which perform well in high dimensional spaces; we’ll discuss one such set of 
techniques, called Markov chain Monte Carlo, later in the course. However, integration 
still tends to be a much more difficult problem than optimization, so if possible we would 
like to formulate our learning algorithms in terms of optimization. Let’s now look at the 
maximum a-posteriori (MAP) approximation, a way of converting the integration problem 
into an optimization problem. 


3.3 Maximum a-posteriori (MAP) approximation 

We worked through two examples of the full Bayesian approach: Bernoulli and Gaussian 
models. In both cases, we saw that as more data is observed, the posterior distribution 
becomes more and more peaked around a single value. This suggests that maybe we can 
get away with summarizing the posterior with a single point estimate. The maximum a- 
posteriori (MAP) approximation chooses the parameters which are most likely under 
the posterior, i.e. 

#map = arg max p(9 \ V) 

0 

= arg max p(Q,V) 

0 

= arg max p(6) p(V \ 6) 

0 

= arg max log p(6) + log p(V \ 0) 

0 

Observe that maximizing log p(T> \ 9) is equivalent to maximum likelihood estimation, so 
the only difference between MAP and ML is the addition of the prior term log p(9). The 
prior is therefore somewhat analogous to a regularize! - . In fact, if p(9) is a Gaussian 
distribution centered at 0, you get L2 regularization! 

Example 5. Let’s return to our coin flip example. The joint probability is given by: 

log p(9, V) = logp($) + log p(V | 9) 

= const + (a — 1) log 9 + (b — 1) log(l — 9) + Nh log 9 + Nt log(l — 9) 

= const + (Nh + a — 1) log 9 + (Nt + 6—1) log(l — 9) 


(Here, const is a shorthand for terms which don’t depend on 9.) Let’s maximize this 
by finding a critical point: 


d 
d 9 


log p(9,V) 


N h + a-l 
9 


Nt + 6 — 1 

1 - 9 


Setting this to zero, we get 
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We can summarize the results of the three different methods in the following table, 
for a = b = 2. 



Formula 

Nh = 2, N-p = 0 

Nh = 55, Nt = 45 

0ML 

N h 

Nh+Nt 

1 

m = °- 55 

0pred 

N H +a 

| « 0.67 

m « °- 548 

0MAP 

N H +a— 1 

| = 0.75 

m * °- 549 

N —2 


When we have 100 observations, all three methods agree quite closely with each other. 
However, with only 2 observations, they are quite different. 0 ml = L which as we 
noted above, is dangerous because it assigns no probability to T, and it will have a 
test log-likelihood of —oo if there is a single T in the test set. The other methods 
smooth the estimates considerably. MAP behaves somewhere in between ML and 
FB; this happens pretty often, as MAP is a sort of compromise between the two 
methods. 

Example 6. Let’s return to our Gaussian example. Let’s maximize the joint prob¬ 
ability: 

1 | w 

log p(n,V) = const - — 2 -(/i - /ipri) 2 - ^ yZ(x M - n) 2 

2 V 2(7 U 

d 1 i N 

— log p{n,V) = - 2-(M - Ppri) + 2 

dp Vi wT 

When we set this to 0, we get exactly the same formula for /Imap as we derived earlier 
for //post- This doesn’t mean the two methods make the same predictions, though. 
The two predictive distributions have the same mean, but the MAP one has stan¬ 
dard deviation <tmap = c, compared with (T pre d = \J u post + cr 2 for the full Bayesian 
approach. In other words, the full Bayesian approach smooths the predictions, while 
MAP does not. Therefore, the full Bayesian approach tends to make more sensible 
predictions in the small data setting. A comparison of the three methods is shown 
in Figure 4. 


3.4 Is MAP a good approximation? 

In both of the examples we looked at, ML, FB, and MAP all made very similar predictions 
in the large data regime, but very different ones in the small data regime. Which setting 
is more typical in practice? 
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Figure 4: Comparison of the predictions made by the ML, FB, and MAP methods about 
future temperatures. (Left) After observing one training case. (Right) After observing 
7 training cases, i.e. one week. 


On one hand, we typically use a lot more data than we did in these toy examples. 
In typical neural net applications, we’d have thousands or millions of training cases. On 
the other hand, we’d also have a lot more parameters: typically thousands or millions. 
Depending on the precise dataset and model architecture, there might or might not be a 
big difference between the methods. 

3.5 Can the full Bayesian method overfit? 

We motivated the Bayesian approach as a way to prevent overfitting. It’s sometimes 
claimed that you can’t overfit if you use the full Bayesian approach. Is this true? In a 
sense, it is. If your prior and likelihood model are both accurate, then Bayesian inference 
will average the predictions over all parameter values that are consistent with the data. 
Either there’s enough data to accurately pinpoint the correct values, or the predictions will 
be averaged over a broad posterior which probably includes values close to the true one. 

However, in the presence of model misspecification, the full Bayesian approach can 
still overfit. This term is unfortunate because it makes it sound like misspecification only 
happens when we do something wrong. But pretty much all the models we use in machine 
learning are vast oversimplifications of reality, so we can’t rely on the theoretical guarantees 
of the Bayesian approach (which rely on the model being correctly specified). We can 
see this in our Toronto temperatures example. Figure 5 shows the posterior predictive 
distribution given the first week of March, as well as a histogram of temperature values for 
the rest of the month. A lot of the temperature values are outside the range predicted by 
the model! There are at least two problems here, both of which result from the erroneous 
i.i.d. assumption: 

• The data are not identically distributed: the observed data are for the start of the 
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Figure 5: The full Bayesian posterior predictive distribution given the temperatures for 
the first week, and a histogram of temperatures for the remainder of the month. Observe 
that the predictions are poor because of model misspecification. 


month, and temperatures may be higher later in the month. 

• The data are not independent: temperatures in subsequent days are correlated, so 
treating each observation as a new independent sample results in a more confident 
posterior distribution than is actually justified. 

Unfortunately, the data are rarely independent in practice, and there are often systematic 
differences between the datasets we train on and the settings where we’ll need to apply the 
learned models in practice. Therefore, overfitting remains a real possibility even with the 
full Bayesian approach. 

4 Summary 

We’ve introduced three different methods for learning probabilistic models: 

• Maximum likelihood (ML), where we choose the parameters which maximize the 
likelihood: 


0ml = argmax£(0) = arg maxlogp(P | 0) 

0 0 

Sometimes we can compute the optimum analytically by setting partial derivatives 
to 0. Otherwise, we need to optimize it using an iterative method such as SGD. 
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• The full Bayesian (FB) approach, where we make predictions using the posterior 
predictive distribution. To do this, we condition on the data and integrate out the 
parameters: 

P {V | v) = f P {d | v)p(V | e) d e. 

Sometimes there’s a closed-form solution to the integral, but otherwise we need to 
solve a difficult high-dimensional integration problem. This is what makes FB so 
difficult to apply in practice. 

• Maximum a-posteriori (MAP), a compromise between ML and FB. We approximate 
the posterior distribution with a single value #map> which maximizes the posterior 
probability 


#map = argmax log p(6 \ V) = argmax logp(0) + log p(V \ 0). 

0 0 

This is similar to ML in that it’s an optimization problem, and the prior term logp(0) 
is analogous to a regularization term. 

All three approaches behave similarly in the setting where there are many more data points 
than parameters. However, in settings where there isn’t enough data to accurately fit the 
parameters, the Bayesian methods have a smoothing effect, which can result in much more 
sensible predictions. Later in this course, we’ll see models where each of these methods is 
useful. 
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